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SN3-91— 19 March 1, 1991 

TO: Distribution 

FROM: SN3/E.L. Christiansen 

SUBJECT; Whipple Shield Sizing Equations 

This memorandum provides Whipple shield sizing equations which 
are recommended for use in current Space Station Freedom (SSF) 

Work Package 2 (WP-2) trade study activities. These equations 
are modifications of the 1969 Cour-Palais predictor equations 
which were used by McDonnell Douglas to size WP-2 meteoroid and 
debris shielding for the SSF preliminary design review. Recent 
hypervelocity impact (HVI) test results have shown that the 
original 1969 Cour-Palais predictor improperly scales to the 
particle sizes that the WP-2 shields must be designed to protect 
against. The original equations have been modified to correct 

this scaling deficiency. Substantial increases to WP-2 shielding 
weight estimates are indicated by the modified equations. 

Several possibilities exist, however, to reduce the weight of 
WP-2 shielding (refer to memo SN3-91-25) . These equations will 
be updated in the future as warranted by the results of ultra- 
high speed (>10 km/sec) HVI tests and further analyses. 

BACKGROUND 

For the preliminary design review (PDR) , McDonnell Douglas Space 
Systems Company (MDSSC) baselined an aluminum Whipple (2-sheet) 
shield for meteoroid/debris protection of all WP-2 critical 
equipment (Ref.: McDonnell Douglas, MDC H4807, Structural 

Integrity Report, June 1990). MDSSC specified a 4.0 inch 
standoff distance and 2219-T62 aluminum for the PDR shields. A 
simplified method was used by MDSSC to size the thicknesses of 
the bumper and rear wall of the Whipple shields. A "design" 
particle size was calculated for each surface of a critical 
element from probability of no-failure requirements, environment 
models, surface area, and orientation considerations. The bumper 
thickness was determined by assuming it was 20% of the "design" 
particle diameter; i.e. , 

t b = 0.2 d (1) 

1969 Cour-Palais Equation 

The 1969 Cour-Palais predictor equation ("non-optimum") was used 
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by MDSSC to derive rear wall thickness of PDR shielding. This 
predictor is given by Equation 2 (Ref. Cour-Palais: "Meteoroid 
Protection by Multiwall Structures," AIAA Paper No. 69-372, 

1969) . 

t M = 0.055 (S p S b ) 1/6 M 1/3 V n /S 0 - 5 (70/a) 0,5 (2) 

Nomenclature 


d projectile diameter (cm) 

S density (g/cc) 

M projectile mass (g) 

S spacing between bumper and rear wall (cm) 

a rear wall yield stress (ksi) 

t thickness (cm) 

© impact angle (deg) measured from surface normal 

V projectile velocity (km/ sec) 

V n normal component of pro j . velocity (km/sec) = V cos 0 

Subscripts: 

b bumper 

p projectile 

w rear wall 


The 1969 Cour-Palais Whipple equation was used in the Apollo 
program to extrapolate test data to meteoroid impact conditions. 
Meteoroid hazard assessments for the Apollo command module and 
lunar module were made using this equation. 

HYPERVELOCITY IMPACT DATA 

The 1969 Cour-Palais equation was derived from HVI data for 
particle diameters up to -0.16 cm. However, it has been shown to 
be unconservative when scaling to particle diameters that the 
WP-2 shields must protect against (-1 cm) . Table 1 contains a 
large historical HVI database for Whipple shields that has been 
updated with recent HVI test results. Projectiles up to 1.9 cm 
(10 g) are in the database. The following summarizes the range 
of HVI data in Table 1: 


Parameter 


Range 


d projectile diameter 
$ projectile density 
$ b bumper density 
V n Velocity 
t b /d Bumper th'k to dia. 
S/d Spacing to dia. 
a Yield Stress 


0.04-1.9 cm 

1.14-2.8 g/cc (nylon, glass, Al) 
2. 7-2. 8 g/cc (aluminum alloys) 
6. 5-7. 5 km/sec (majority) 
0.08-0.64 
13-96 

18-70 ksi (aluminum alloys) 
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The Table 1 data set has been selected to be at or near the 
ballistic limit (i. e. , perforation or spall threshold) of the 
rear wall. Of the 55 data points in Table 1, 23 are on the 
failure side of the ballistic limit (i.e., rear walls with 
perforations or detached spall) and 32 have thick enough rear 
walls to prevent failure. 

Table 2 compares predictions made by the original Cour-Palais 
equation (Eq. 2) and the actual experimental rear wall thickness 
for the data set in Table 1, Table 3 summarizes Table 2 results 
and illustrates that the original 1969 Cour-Palais equation 
predicted only 30% of the failures in the data set. None of the 
successful predictions of failure occured with particle sizes 
greater than 0.2 cm. This implies that application of the 1969 
Cour-Palais equation to particle threats of greater than 0.2 cm 
diameter, as was done to develop the WP-2 PDR shielding designs, 
will result in underestimating the required HVI shielding. 


Table 3. 1969 Cour-Palais Equation Predictions 

1969 Cour-Palais 


* 

Predicted Failures 7 

Prediction Accuracy 30.4% 

Mean tw calc/tw exp't 0.82 

Standard Dev. tw calc/tw exp't 0.81 

Predicted No-Failures 27 

Prediction Accuracy 84.4% 

Mean tw calc/tw exp't 0.63 

Standard Dev. tw calc/tw exp't 0.55 


* Predicted failures out of the failures in Table 1: 23 data points 
** Predicted no-failures out of the no-failures in Table 1: 32 data points 


MODIFIED COUR-PALAIS WHIPPLE EQUATIONS 

The HVI data in Table 1 has been used to modify the original 
Cour-Palais equations as given below (equations 3 and 4). 


tb 

= 

0.25 d S p /S b 

(3) 

tw 

= 

C d 0,5 (S p S b ) 1/6 M 1/3 V n /S 0,5 (70/a) 0-5 

(4) 

* • 2 2/3 

Where coefficient C = 0.16 cm -sec/g ' -km. 



Bumper thickness, in Equation 3, has been adjusted by the ratio 
of projectile to bumper density. This change will result in 
reductions in the MDSSC weight estimate for bumpers on surfaces 
of critical equipment that are only exposed to the meteoroid flux 
(because meteoroids are low density) . The coefficient in 


3 





Equation 3 has been increased to 0.25 from the MDSSC PDR approach 
of using 0.20. This is required to reduce the possibility of 
underestimating the required rear wall thickness with small 
standoff distances (i.e. , when S/d < 15). If standoff distance 
is increased (S/d > 30), then the original 0.20 coefficient can 
be substituted without reduction in accuracy. 

Equation 4 is a slightly modified version of the 1969 Cour-Palais 
equation. The coefficient in the original equation (0.055) has 
been replaced with a diameter scaled coefficient (0 . 16*d 0-5 ) to 
reflect results of data from large particle hypervelocity impact 
testing. These equations are derived from the analysis of Cour- 
Palais. Data from Table 1 was used to derive the diameter-scaled 
coefficient in Equation 4. The "K" factor (K = C * d 0,5 ) curve 
fit to the data is shown in Figure l. 

As given in Table 4, the modified Cour-Palais equation more 
accurately predicts failure; successfully predicting over 65% of 
the failure cases. All 8 data points that Equation 4 did not 
successfully predict to fail were very marginal failures 
(detached spall without perforation or small perforations) ; and 
only 2 of the 8 are for larger particles (0.5 cm and larger). On 
the other hand, the original 1969 Cour-Palais equation 
successfully predicts only 7 of 23 failures; and all of these are 
for small particles (0.2 cm and smaller). 


Table 4. Modified Cour-Palais Equation Predictions 

Modified COur- Palaia 


Predicted Failures 


15 

Prediction Accuracy 


65.2% 

Mean tw calc/tw exp* t 


1.15 

Standard Dev. tw calc/tw 

exp 1 1 

0.45 

** 

Predicted No-Failures 


23 

Prediction Accuracy 


71.9% 

Mean tw calc/tw exp* t 


0.90 

Standard Dev. tw calc/tw 

exp • t 

0.36 


* Predicted failures out of the failures in Table 1: 23 data points 
** Predicted no-failures out of the no-failures in Table 1: 32 data points 


A more conservative 0. 19*d°‘ 5 coefficient was originally proposed 
(Memo SN3-90-13 1 , Christiansen: "Shield Sizing Equations", 
10/12/90) to guarantee that the rear wall would not fail with 
V r) > 6 . 5 km/sec. However, projectile fragmentation and melting 
will increase between 6.5 and 7 km/sec. This will result in less 
rear wall damage as velocity increases from 6.5 to 7 km/sec. 

When "K" factors from HVI data shots with velocities less than 7 
km/sec were adjusted to 7 km/sec, the proposed 0.16*d 0-5 
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coefficient better fits the data (Figure 2) . 

Constraints; 


The ballistic limit calculated by the modified equations is a 
local deformation without detached spall or fracture; that is, 
the no-failure criterion is no perforation or detached spall of 
the rear wall (corresponding to damage categories: D1-D2, E1-E2, 

F1-F3 as given in Figures 3-5) . Equation 4 will potentially 
under predict the required rear wall thickness to prevent failure 
for normal component velocity (V n ) of less than 7 km/sec, S/d 
ratios of less than 15, and t b /d ratios of less than 0.15. A set 
of ballistic limit equations for aluminum Whipple shields is 
given in memo SN3-91-21 that can be used over the entire debris 
velocity range from 0 to 15 km/ sec. The modified Cour-Palais 
equation presented in this memo is included in the set of 
ballistic limit equations given in SN3-91-21. 

Technical Basis : 

Technical rationale for using the proposed equation include: 

• The approach is a conservative method for predicting the HVI 
response of Whipple shields to impact conditions beyond 
testing capability. Without test data and/or detailed 
analyses, a conservative approach is more appropriate for 
designing the protection systems of critical equipment since 
they are directly associated with crew safety and SSF 
survi vabi 1 ity . 

• The equation is a semi-analytical/empirical approach. It 
uses test data to "anchor" the prediction at the highest 
impact velocities attainable in the laboratory, and uses a 
conservative analytical approach to extrapolate to higher 
velocities. 

/ 

• The basis of the proposed scaling equation is that the 
debris cloud impacting the rear wall will contain solid 
particulates under certain impact conditions (even at 
velocities ranging from 7-15 km/sec) . These solid particles 
are the determining factor in sizing the rear wall. Late 
time fragments from the bumper, impacts of irregularly 
shaped debris particles (such as plates, rods, hollow 
objects, etc) and impacts at highly oblique angles or with 
high t b /d ratios can potentially release solid fragments; 
and are a likely occurrence in the real environment. 

BUMPER Application ; 

In probability analysis tools, such as the modified BUMPER 
program, an equation is needed that relates critical particle 
diameter to fail a given structure with impact velocity and 
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angle. Rearranging Eq.4 results in the following penetration 
equation for the critical particle diameter causing failure: 

d crit = 3.918 t w 2/3 5 p ' 1/3 <s b ' 1/9 V' 2/3 cos' 2/3 0 S 1/3 (a/70) 1/3 (5) 

This equation indicates the ballistic limit of a given structure 
(i.e., the impact conditions causing rear wall failure) scales 
with constant impactor kinetic energy (i.e., 

ballistic limit a K.E. 1/3 a d cHt V n 2 ' 3 S ^ m a constant for a 
given structure) . This is consistent with NASA practice for 
conservatively extrapolating beyond HVI test conditions (Ref. 

NASA SP— 8042 , NASA Space Vehicle Design Criteria: "Meteoroid 
Damage Assessment," May 1970, p. 37) , as well as the presumption 
that the debris cloud will contain solid particulates in a 
significant fraction of real encounters occuring at very high 
impact velocities (> 8 km/sec) and cratering theory for the 
penetration of solid particles into the rear wall. 

An example of the critical particle size calculated by Eq. 5 that 
fails the MDSSC PDR design for the orbital debris shields on the 
propulsion module is shown as a function of velocity in Figure 6 
(all materials are A1 2219-T62) . This figure shows the bumper 
thickness to critical particle diameter ratio ranges from 0.28 to 
0.47. A similar example is shown in Figure 7 for the debris 
shielding on the habitation and laboratory modules. 

DESIGN IMPLICATIONS 

Changing from the original 1969 Cour-Palais equation to the 
modified Cour-Palais equation will more than double rear wall 
thickness and weight for the debris particle sizes that are 
currently driving the design of WP-2 protection systems. Since 
the debris shields are -70% of total shielding weight, an 
increase of approximately 80% to 90% in WP-2 shielding can be 
expected for using the proposed shield sizing equations to meet 
current requirements and environment — given no changes in design 
approach. 

Weight Reduction Strategies 

There are several options that could reduce shielding weight, 
without reducing protection performance. These options are 
presented in another memo SN3-91-25: "Weight Reduction Options 

for Meteoroid/Debris Shielding" . 

COMPARISON WITH WILKINSON PREDICTOR 

The Wilkinson predictor is used by the WP-1 contractor (Boeing) 
to predict Whipple shield performance beyond test conditions. 

The Wilkinson predictor is given by the following equations (Ref. 
J.P.D. Wilkinson: "A Penetration Criterion for Double-Walled 

Structures Subject to Meteoroid Impact," AIAA Journal, Vol.7, 
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No. 10, pp. 1937-1943 , October 1969). 

For [t b «S b /( 5 P d ) 3 > 1 : 

t u = M V/ ( 1 . 44 L 2 S 2 SJ (6) 


For [t b S b /(S p d)] < 1 : 

t H = M 473 5 b 2/3 V/ [1.44 (7T/6) 173 L 2 t b 5 b S 2 SJ (7) 

L, is a material constant (0.425 for Al 2024-T3 , A1 2024-T4, and 
A1 2024-T351; 0.292 for Al 6061-T6; 0.345 for Al 7075-T6; 0.297 
for Al 2014-T6; 0.28 for Al 2219T87; 0.20 for Al 3003-H12, Al 
3 003-H14 , and Al 1100-H14) . 

The Wilkinson equations were applied to predict rear wall 
thickness using the data in Table 1. Table 5 shows that the 
Wilkinson equations are very unconservative for this data set. 

The calculated rear wall thicknesses required to prevent failure 
using Wilkinson were smaller than the actual experimental rear 
wall thicknesses for all 55 HVI tests; indicating no failure is 
expected to occur in any of the 55 cases. Thus, Wilkinson's 
accuracy in predicting shield failure is 0% for this data set. 
Since Wilkinson always predicted shield success, it was 100% 
accurate in predicting no-failure. However, for the no-failure 
data, the rear wall thickness calculated by Wilkinson averaged 
only 8% of the rear wall thickness used in the experiments. This 
is quite optimistic considering the data in Table 1 was 
intentionally selected to be close to the failure point of the 
shield. 

The poor comparison between Wilkinson prediction and experiment 
is not unexpected. The Wilkinson equation was formulated almost 
entirely from analytical considerations, without use of 
experimental data to determine coefficients or material 
constants, and little attempt was made to compare predictions 
with experimental data. 

COMPARISON WITH NYSMITH EQUATION 

For some configurations, WP-4 contractors have used the Nysmith 
equation to extrapolate Whipple shield performance beyond test 
conditions. The Nysmith equation is given by the following 
equation (Refs: C.R. Nysmith: "An Experimental Impact 

Investigation of Aluminum Double-Sheet Structures," AIAA Paper 
No. 69-375, AIAA Hypervelocity Impact Conference, April 30 - May 
2, 1969. C.R. Nysmith: "Pentration Resistance of Double-Sheet 

Structures At Velocities to 8.8 km/ sec, " NASA TN D-4568, May 
1968. ) . 

t w = 5.08 d v 0 - 2778 (t b /d) ' 0,5278 (S/d)* 1 - 3889 ( 8) 
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Table 5 shows the results of using Equation 8 to predict rear 
wall thickness for the data in Table 1. As indicated, Nysmith 
only predicts 13% of the test cases where the shield failed; an 
unconservative result. This is somewhat surprising. The Nysmith 
predictor is purely empirical (HVI test data is used in its 
formulation) . The poor comparison here can be explained by 
noting that Nysmith formulated his equation using HVI data for 
glass projectiles on aluminum targets. Because the density of 
glass is less than aluminum, this correlation would underestimate 
the effect of aluminum projectiles. However, of more concern is 
use of the Nysmith equation to extrapolate beyond test 
conditions. The Nysmith equation is derived strictly from HVI 
data, and contains no theoretical justification for extrapolating 
beyond test conditions. 


Table 5. Comparison of Wilkinson and Nysmith Equations 


* 

Predicted Failures 
Prediction Accuracy 
Mean tw calc/tw exp't 
Standard Dev. tw calc/tw exp't 

Predicted No-Failures 
Prediction Accuracy 
Mean tw calc/tw exp't 
Standard Dev. tw calc/tw exp't 


Wilkinson 

Nvsmith 

0 

3 

0% 

13.0% 

■W 

r-t 

• 

o 

0.58 

0 . 18 

0.39 

32 

32 

100% 

100% 

0.08 

0.37 

0.05 

0.16 


* Predicted failures out of the failures in Table 1: 23 data points 
** Predicted no-failures out of the no-failures in Table 1: 32 data points 


Distribution 

ET13/Ray Nieder MDSSC/Burton G. Cour-Palais 

KC2/Dale Haines 

SN/Doug Blanchard 

SN/ David Thompson 

SN3 /Drew Potter 

SN3 /Jeanne Lee Crews 
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FIGURE 1 

Whipple Shield "K" Factor 

Modified Cour— Palais Whipple Non— Opt Eq 



ja^aDj >| 


Proj. Dia. I 
(U)=target failed 



FIGURE 2 

Whipple Equation "K" Factor 

K=*tw (S*Y/70)~ .5*(pb*p)^— .5*m ~-.33*Vn ~ - 1 
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Class 

CATEGORY A: SINGLE CRATER PATTERN - LOW VELOCITY 

PROJECTILE REMAINS INTACT 

A1 • NO PERFORATION OR REAR SURFACE DEFORMATION 
. CRATER DIAMETER APPROXIMATE SIZE OF PROJECTILE 

A2 • NO PENETRATION 

• CRACKS OR SPLITTING MAY BE PRESENT 
. REAR SURFACE DEFORMATION 


| IMPACT 
I DIRECTION 




A3 • PENETRATION 

• HOLE DIAMETER APPROXIMATE SIZE OF PROJECTILE 


E 77 77 !% 


B1 


CATEGORY B: SINGLE CRATER PATTERN - HYPERVELOCITY 

PROJECTILE REMAINS INTACT 


• NO PERFORATION OR REAR SPALL 

• SINGLE ROUNDED CRATER 

. FRONT SURFACE LIP OR SPALLATION 



B2 • NO PERFORATION, BUT WITH ATTACHED REAR SPALL 

• SINGLE ROUNDED CRATER 

• FRONT SURFACE LIP OR SPALLATION 



B3 • NO PERFORATION, BUT WITH DETACHED REAR SPALL 

• SINGLE ROUNDED CRATER 

• FRONT SURFACE LIP OR SPALLATION 

• LIGHTTIGHT 



B4 • PERFORATION DUE TO CRATER AND REAR SPALL MEETING 
(HOLE DIAMETER < 2 mm) 

• FRONT SURFACE LIP OR SPALLATION 
•NOT LIGHT TIGHT 




B5 •PENETRATION 

• HOLE FORMED BY CRATER AND DETACHED SPALL 
(HOLE DIAMETER > 2 mm) 

• FRONT AND REAR SURFACE LIPS OR SPALLATION 


VZ2Z/& 


figure 3. Damage Classification for Shielded Metallic Targets 


REF. Dahl and Cour-Palais* "stanrisrHi,^' 

°“ a ?f. C1 ? ssifica Wcn and ImP 

Metallic Targets " , 1990 xor 



Class 

CATEGORY C: MULTIPLE CRATER PATTERN 
PROJECTILE BREAKS UP INTO COARSE AND FINE FRAGMENTS 

Cl • NO PERFORATION, REAR SURFACE DEFORMATION OR SPALL 

• RANDOM SURFACE CRATERS, PITTING OR EROSION 

C2 • NO PERFORATION, BUT WITH ATTACHED SPALL(S) OR REAR 
SURFACE DEFORMATION 

• RANDOM SURFACE CRATERS, PITTING, OR EROSION 

C3 • NO PERFORATION, BUT WITH DETACHED SPALL(S) 

• RANDOM SURFACE CRATERS, PITTING OR EROSION 
•LIGHTTIGHT 

C4 • PERFORATION 

• CRACKS OR SMALL HOLE(S) (ALL HOLE DIAMETERS < 2 mm) 

• NOT LIGHTTIGHT 

C5 • PENETRATION 

• LARGE HOLE(S) (APPLICABLE IF ANY HOLE DIAMETER > 2 mm) 


I IMPACT 
I DIRECTION 







CATEGORY D: CENTRAL CRATER PATTERN 

PROJECTILE BREAKS UP INTO FINE PARTICLES 

D1 • NO PERFORATION OR REAR SPALL 

• CENTRAL SURFACE CRATER, PITTING OR EROSION 

D 2 • NO PERFORATION, BUT WITH ATTACHED SPALL 

• CENTRAL SURFACE CRATER, PITTING, OR EROSION 

D3 • NO PERFORATION, BUT WITH DETACHED SPALL 

• CENTRAL SURFACE CRATER, PITTING OR EROSION 
•LIGHTTIGHT 

D4 • PERFORATION 

• CRACKS OR SMALL HOLE(S) DUE TO CRATER AND SPALL 
MEETING (ALL HOLE DIAMETERS < 2 mm) 

• NOT LIGHT TIGHT 

D5 • PENETRATION 

• LARGE HOLE(S) FORMED BY CRATER AND DETACHED SPALL 
(APPLICABLE IF ANY HOLE DIAMETER > 2 mm) 







figure 4. Damage Classification for Shielded Metallic Targets (Cont.) 

x 

REF. Dahl and Cour-Palais: "Standardization of Impact 

Damage Classification and Measurements for 
Metallic Targets", 1990 



Class 


IMPACT 

DIRECTION 


CATEGORY E: RING CRATER PATTERN 
PROJECTILE BREAKS UP INTO VERY FINE PARTICLES f 


El • NO PERFORATION OR REAR SPALL 

• RING CRATERS SURROUND CENTRAL SURFACE CRATER, 
PITTING, OR EROSION 

E2 • NO PERFORATION 

• RING CRATERS WITH SPALL PIMPLES ATTACHED AND/OR 
CENTRAL SPALL ATTACHED 

• CENTRAL SURFACE CRATER, PITTING, OR EROSION 

E3 • NO PERFORATION 

. RING CRATERS WITH SPALL PIMPLES DETACHED AND/OR 
CENTRAL SPALL DETACHED 

• CENTRAL SURFACE CRATER, PITTING, OR EROSION 
•LIGHTTIGHT 





E4 • PERFORATION _ _ _ « 

• HOLE(S) DUE TO CRATER(S) AND SPALL(S) MEETING W2}QjF sSs>£ ^Zh ^ZZA 

•NOT LIGHT TIGHT 


E5 * PENETRATION 

• LARGE HOLE PUNCHED OUT DUE TO RING PERFORATIONS *773 tZZS- 

AND IMPULSIVE LOAD 


CATEGORY F: NON-PARTICULATE IMPULSIVE LOADING 

PROJECTILE BECOMES MOLTEN LIQUID OR VAPOR 

FI • NO PERFORATION OR REAR SPALL 

• SURFACE PITTING OR MOLTEN SPLASH 

F2 • NO PERFORATION 

• SPALL PRESENT, ATTACHED OR DETACHED 
. SURFACE PITTING OR MOLTEN SPLASH 

F3 • NO PERFORATION 

• DENTED, BUT INTACT 

• SURFACE PITTING OR MOLTEN SPLASH 

• LIGHTTIGHT 

F4 • PERFORATION 

• DENTED AND SPLIT 

• SURFACE PITTING OR MOLTEN SPLASH 

• NOT LIGHT TIGHT 

F5 • PENETRATION BY IMPULSIVE LOAD FAILURE 

• PETALLED HOLE 

• SURFACE PITTING OR MOLTEN SPLASH 


'/ 




figure 5 . Damage Classification for Shieided Metallic Targets (Cont.) 


REF. Dahl and Cour-Palais: "Standardization of Impact 

Damage Classification and Measurements for 
Metallic Targets", 1990 



Critical Particle Size vs. Velocity 

MDSSC PDR, .052" bump, S=4", .099" wall 
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